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ABSTRACT 

Optical sensors aboard Earth orbiting satellites such as the next generation Visible/Infrared Imager/Radiometer Suite 
(VIIRS) assume that the sensors’ radiometric response in the Reflective Solar Bands (RSB) is described by a 
quadratic polynomial, in relating the aperture spectral radiance to the sensor Digital Number (DN) readout. For 
VIIRS Flight Unit 1, the coefficients are to be determined before launch by an attenuation method, although the 
linear coefficient will be further determined on-orbit through observing the Solar Diffuser. In determining the 
quadratic polynomial coefficients by the attenuation method, a Maximum Likelihood approach is applied in carrying 
out the least-squares procedure. Crucial to the Maximum Likelihood least-squares procedure is the computation of 
the weight. The weight not only has a contribution from the noise of the sensor’s digital count, with an important 
contribution from digitization error, but also is affected heavily by the mathematical expression used to predict the 
value of the dependent variable, because both the independent and the dependent variables contain random noise. In 
addition, model errors have a major impact on the uncertainties of the coefficients. The Maximum Likelihood 
approach demonstrates the inadequacy of the attenuation method model with a quadratic polynomial for the 
retrieved spectral radiance. We show that using the inadequate model dramatically increases the uncertainties of the 
coefficients. We compute the coefficient values and their uncertainties, considering both measurement and model 
errors. 
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1. INTRODUCTION 

As the next generation remote sensing optical instrument to the MODerate Resolution Imaging 
Spectroradiometer 1 (MODIS), a VIIRS instrument 2 will be carried aboard each platform of the Joint Polar Satellite 
System and the National Polar-orbiting Environmental Satellite System Preparatory Project (NPP). The VIIRS 
instrument offers high quality imaging capabilities in visible and infrared bandwidths with slightly finer spatial 
resolution than the MODIS’, providing more accurate global weather and environmental data. Due to limited on- 
orbit calibration means, to accurately measure the detected spectral radiance, extensive prelaunch efforts have been 
conducted to calibrate and characterize the sensors. 

To facilitate aperture spectral radiance determination, we assume that the radiance relates a sensor’s DN count 
through a mathematical expression. In order to gain robustness and accuracy, extending the linear mathematical 
form used in the MODIS for RSB and taking the similar expression used for both the MODIS and VIIRS Thermal 
Emissive Bands 3,4 , we use a quadratic polynomial for the radiance 



Z det = c 0 + c \dn + c 2 dn 2 , 


( 1 ) 


where dn is the background subtracted Digital Number. The polynomial coefficients are determined pre-launch, 
using a National Institute of Standards and Technology (NIST) traceable Spherical Integrating Source (SIS-100). 
Because the SIS-100 light source is not stable enough over time and its output radiance is not accurately known with 
a relative uncertainty of a few percent, an attenuation method is applied to determine the c-ratios: c 0 /c, and c 2 / c, . 

In the attenuation method, the sensor measures the same output radiance of SIS-100 with and without an 
attenuator . Since the time duration between reading the digital counts is short (2-3 minutes), it is reasonable to 
assume that the SIS-100 output radiance does not change much in the time duration. Because the attenuator is an 
opaque plate with small holes to allow light through, the attenuator’s transmittance is unchanged over the entire SIS- 
100 output radiance levels. Consequently, the ratio of the detected spectral radiances with and without the attenuator 
may be described by 

_ c 0 /c l +dn in +(c 2 /c l )dnl 

^ ( \ 2 9 ^ 

c o / c i + dn out + \c 2 / c l )dn out 

where dn m and dn out denote the dn with and without the attenuator, respectively, and r the transmittance, a 
constant across lamp levels. The c-ratios are determined through Eq. (2). 

The legacy approach to determine the c-ratios uses a least-squares process with equal weights, namely, each data 
point in the least-squares has the same weight. This approach can be problematic, especially when most of the 
measurement data points are poorly determined, resulting in inaccurate values for the parameters. In addition, the 
equal-weight approach does not give information on whether the model Eq. (2) is adequate. model is correct. An 
inadequate model, if forced to be used, will have a significant impact on the fitting parameters. 

To address the problems of the legacy approach, in this study, we use a Maximum Likelihood 5 approach. The 
Maximum Likelihood assigns appropriately more weight to those data points with smaller measurement 
uncertainties so that with the determined c-ratios the likelihood of observing the data is at its maximum. Unlike the 
equal-weight approach, the Maximum Likelihood approach gives us a condition to determine the validity of the 
model with a very high probability. If the model is inadequate, we shall calculate the impact on the c-ratios due to 
the inadequacy. 

In the next section, we provide our theoretical description of using the Maximum Likelihood to determine the c- 
ratios. We describe a method for computing c-ratio errors due to model inadequacy. In Section 3, we give an 
example of the results obtained by applying the theories developed in Section 2 to VIIRS FU 1 Thermal Vacuum 
(TV) data for RSB. We use the results to determine whether the quadratic polynomial model Eq. (2) adequately 
describes our data. In Section 4, we try to gain insight for the procedure given in Section 2. Specifically, we consider 
the impact of model inadequacy on fitting parameter uncertainties. Finally, in Section 5 we summarize the results of 
this study. 

2. MAXIMUM LIKELIHOOD APPROACH 

To use existing data regression functions, we need to write Eq. (2) into a form of y = f(x; p) , where all the 
fitting parameters are represented by p. Let x = dn out and y = dn m , Eq. (2) can be written as 



( 3 ) 



where z = c 0 / q + x + (c 2 / Cj)x ? . The background DN is averaged over 48 consecutive Half Angle Mirror (HAM) 
angular positions (denoted as the sample positions) over the Space View port for each HAM side and each scan, and 
the Earth View DN is the average over 25 scans with repeated 3-sigma outlier rejections at each of the 120 angular 
positions (since a 3 -sigma rejection changes the mean as well as its standard deviation, it is necessary to run the 3- 
sigma rejection repeatedly until the final results stabilize). Due to the DN digitization, a zero standard deviation of 
the mean is possible. When a zero standard deviation occurs, we assign a quantization error as the standard deviation 
of the averaged DN. 

Even if the model is accurate, Eq. (3) may not hold due to measurement errors in y and/or x. The errors in y 
and/or .r make the Maximum Likelihood approach viable and the weight is the inverse of the variance of 
y - /(. x; p) where f(x; p ) stands for the mathematical expression on the right hand side of Eq. (3) and p denotes 
the c-ratios. 

In order to compute the Maximum Likelihood weight, we assume that the measurement errors follow Gaussian 
distributions. As a result, the probability distribution for y-f(x; p) is also Gaussian. The variance of y - f(x; p) 
is computed by 

cr 2 = var(y)+ var (f(x; /?))-2covar(>>, f(x\ p)) . (4) 

Since the errors in dn [n and dn out are not correlated and f(x; p)& rxdn out , Eq. (4) is simplified to become 

a 2 = a 2 (dn in )+ T 2 a 2 (dn out ) . (5) 

Eq. (4) indicates that the weight in the Maximum Likelihood approach depends on the mathematical form of the 
fitting function. 

In order to avoid any adverse effects on the determined c-ratios from the variation of the attenuator’s 
transmittance over sample positions, we carry out our fitting procedure at each sample position. The detailed 
procedure is described in the following steps. 

(a) At each sample position, the measured data points are denoted as ( dn m , dn out ) sample , for all valid SIS- 100 
radiance levels. Using a data regression procedure, we find the parameters (c 0 /c l ,c 2 /c l ) sample which minimize 


X [dn m (n)- f(dn out («); p)f l a 2 {n) , 


where N is the number of SIS- 100 radiance levels. The data regression procedure we use is the Interactive Data 
Language function MPFIT( . . .) 6 which performs a Levenberg-Marquardt least-squares procedure with additional 
features such as allowing setting limits on the values of fitting parameters. We use repeated 4-sigma rejections until 
the rejection stabilizes. 



(b) We determine the simple means of (c 0 / q, c 2 / q) sample over the sample positions with repeated 3-sigma outlier 
rejections for a set of 120 2-element vectors. At this stage, the means of (c 0 / q , c 2 / q ) sample are for a single HAM 
side. 

(c) Because the c-ratios in VIIRS RSB are not HAM side dependent (assuming the mirror does not contribute to 
inter-detector cross talks), we compute the simple mean of ^c 0 /q, c 2 / c i) sample over HAM sides to obtain 

(c 0 /q,c 2 /q) sample>HAM . 


At the end of step (c), we must address the question of whether the quadratic model is adequate. If the model is 
adequate, we can determine the covariance matrix for the c-ratios in step (b) (from the 120 c-ratio covariance 
matrixes). If the model is inadequate, indicated by small goodness-of-fit numbers, we have three choices to improve 
our capability to determine the aperture radiance. One is to change the model. Another is to divide the fitting region 
into two or more smaller ones and carry out the fitting process separately in those regions. The third is to apply the 
inadequate quadratic polynomial for the entire region, in order not to increase algorithm complexity. 

Using an inadequate model can generate large differences over some radiance regions between the measured and 
the fit-obtained values in r . The large differences in r can be quantified by the variance of r . This variance 
enlarges the uncertainties of (c 0 / q , c 2 / c \ ) sam pie ham • Details to compute these uncertainties are described in step 

(d). 

(d) We denote the simple mean of the r differences over samples and HAM sides at a particular radiance level as 
l At) , TT>W , . (At) , TT>W , is the sum of the differences induced from random noise and model 
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Since Ar modelrad and ( Ar ) sampledlAMirad noise 
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are statistically independent, Eq. (6) gives 
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In Eq. (7), z; is a positive number, equal to 1 if the digital counts from consecutive samples from the same scan are 
not correlated. To determine ij , we can average the digital counts over the samples first and then take the scan 
average. In Eq. (8), on the left hand side, the outer <...> indicates the statistical average and 7/ samp i e = 120 . 

To compute the variance in r over the entire radiance range, we realize that the contributions to the variance 
from the random noise and the model errors depend differently on the number of radiance levels. The contribution 
from the random noise is 

var noise ( T ) “ D 'Z.l f } samplej HAM radj noise j 2 ^ > (9) 



which is roughly proportionally to 1 IN, indicating that as long as we have enough number of radiance levels, the 
contribution to the variance of r can be very small. The model error generated r difference, however, contributes to 
the t variance in a mean squared fashion. Therefore, when model errors dominate, we can compute the overall r 
variance as 

var(r) = var model (r) + var noise (r) , (10) 

where the model error created r variance var model (r) is defined as 

Lmax 

var model (r)= Jvar model (r, L)dL/(L ma3i -L min ) . (11) 

Anin 


In Eq. (11), the model error created variance of r at radiance L is 


var model (r, L) 
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Fig. 1 . Each data point is from a fit of Eq. (3) to measurement data at a fixed attenuator transmittance, (a) c 0 / q versus 
T . (b) (?2 / c \ versus T . Results are for detector 1 of band Ml at (F1AM side A, Fixed High Gain, Nominal temperature 
plateau, control electronics side B). 

To use the variance of r to compute the uncertainties of (c 0 / c t , c 2 / c, ) , we rely on the findings between the 
changes in ( c 0 1 c x , c 2 I c x ) and the change in r . Mathematically, the relations can be visualized through Eq. (2) 
which gives that the impact of a small change in r on the c-ratios obeys the following relation 
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(13) 





Eq. (13) demonstrates that for a small change in r at a particular radiance level, the changes in the c-ratios fall on a 
straight line. Taking account all the radiance levels, in a minimum error sense, such as least-squares, we find the 
changes in the c-ratios due to a small change in r behave as S(c 0 / cq ) oc (dn out ^ x St and S(c 2 / c { )cc -St /{cln 0ut 'j . 

To visually show the relation between the change in r and the changes in the c-ratios, we use Eq. (3) and carry out 
data regressions at a few fixed r at sample position one for detector 1 of band Ml (412 nm) at (ElAM side A, Fixed 
High Gain, Nominal temperature plateau, control electronics side B). Figs, la and lb show plots for c 0 / c 1 versus r 
and c 2 / C| versus r , respectively. The figures clearly indicate linear relationship between the c-ratios and r . As a 
result, we compute the uncertainties of the c-ratios by multiplying the square root of var(r) by the absolute values 
of the slopes exhibited by the data points in Figs, la and lb, respectively. Furthermore, the figures reveal that if 
model error dominates over the random measurement errors, the correlation coefficient for the c-ratios is -1. 

That we are interested in the variances of the c-ratios is because we want to determine the relative uncertainty of 

c 0 / q + dn + c 2 I C\ x dn 2 , computed by 
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3. RESULTS 

We carry out the steps mentioned in Section 2 and obtain the results for all detectors in RSB, at electronics gain 
settings of Fixed Low and High Gains, for temperature plateaus of Cold, Nominal and Hot, and at control electronics 
sides A and B, using VI1RS TV RC-02 data. The values of the c-ratios agree with those obtained previously with the 
equal-weight approach within the standard deviations from the equal-weight approach. 

To determine the c-ratio variances and the correlation coefficient, after we carry out the data regression at the per 
sample stage, we compute the chi-square and the chi-square at goodness-of-fit of 0.001. The chi-squares for detector 
1 of band Ml at (Fixed High Gain, Nominal temperature plateau, control electronics side B) are about 250 for all the 
samples and HAM sides, significant larger than 145 which is about the chi-square at goodness-of-fit of 0.001, 
indicating strongly that the quadratic model is inadequate. We therefore use step (d) to compute the variances of the 
c-ratios through the variance of r and assign -1 to the c-ratio correlation coefficient. 

In Figs. 2a and 2b, we show our results in diamonds for the c-ratios and their standard deviations for the 
detectors in band Ml at (HAM side A, Fixed High Gain, Nominal temperature plateau, control electronics side B). 
In the figures, error bars stand for the standard deviations. The squares are for the equal-weight approach without 
considering model errors. For the Maximum Likelihood approach without considering model errors, namely without 
step (d), the standard deviations for the c-ratios are about an order of magnitude smaller than the ones from the 
equal-weight approach. After considering model errors, the standard deviations for the c-ratios increase dramatically 
as shown in the figures. 
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Fig. 2. Diamonds are for HAM and sample averaged (simple average) c-ratios from the Maximum Likelihood fits of 
Eq. (3) to measurement data for all detectors in band Ml at Fixed High Gain and Nominal temperature plateau at 
control electronics side B. Bars stand for the standard deviations. Squares are for the equal-weight approach, (a) Cg/q . 

(b) c 2 /q. 



Fig. 3. Relative radiance errors, in percentage, computed from Eq. (13) for band Ml detector 1, at (HAM side A, Fixed 
High Gain, Nominal temperature plateau, control electronics side B). 

With the variances and the covariance of the c-ratios, we compute the relative radiance errors, determined by Eq. 
(14). We show the results in Fig. 3 for detector 8 in band Ml at (HAM side A, Fixed High Gain, Nominal 
temperature plateau, control electronics side B). The dashed vertical lines indicate the dn at the requirement 






specified minimum, typical, and maximum radiances. From the figure, we can see that the relative error is quite 
small (less than 0.3%) over [Z min , Z max ] . When dn becomes smaller than 300 (<dnMin)), with decreasing dn, the 
relative error becomes larger quickly. 


4. DISCUSSION 

To justify the Maximum Likelihood approach and the application of Eq. (10) to compute the variance of t in 
the case of model inadequacy, we consider fitting the model of a constant. We assume that at {x m ; m = 1, • • • , M } , the 
measured value set is {y m }. With a measurement error dy m and a bias function b{x), we have the following 
equation relating the measurement value with the true value 


y m =C + b(x m )+5y m . 

The Maximum Likelihood approach yields that 


(15) 


M 


I \}’m -HO] /( 


c = 


n = 1 


M 

YMa] 

m = 1 


(16) 


We force a constant model of 
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by removing the bias function to obtain 
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Comparing Eqs. (16) and (18), we know 
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Given Eqs. (15-19), the question we need to answer is how we can use the fit-obtained C’ to predict the true 
C + Z>(x) with some level of certainty. 

M M 

We can see from Eq. (19) that C and C’ differ by the weighted sum of the biases 'E J b(x m )/ a m / XU a m . 

m = 1 m = 1 

Although we do not exactly know the bias function, it can be computed approximately. If model errors dominate 
measurement errors, indicated by a very small goodness-of-fit number, we can replace b(x m ) by y m - C' and obtain 



( 20 ) 
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Therefore, C and C’ essentially are the same. 

Next, we need to find the standard deviation of C’ in the presence ofbiases. Assuming covar(y ; ., y •; i ^ /)= 0 , 
the variances of C ’ from the random noise is computed by 

var n 0 ise (C) = -Jf- , (21) 

si lal 

m = 1 


indicating that the variance of C ’ from measurement errors can be very small, as long as we have enough number of 
data points. Consequently, Eq. (21) can not address the effect of the biases since the biases are not random and more 
data points should not decrease the variance of C’. Asa result, in the presence ofbiases, obtaining the variance from 
the output covariance matrix generated by data regression routines will not be correct, whether to use Maximum 
Likelihood weights or a uniform weight. With biases, the variance of C ’ may be computed by 

var(C') = std 2 (b)+ var noise (C' ) . (22) 

std(b) is defined as that over [x min , x max ] , 68.2% of |/i(x)| are smaller than std(b). Eq. (22) may be approximated 
by 

var(C') = J b 2 (x)dx /(x max - x min )+ var noise (C') . (23) 

■^min 

Since we do not exactly know the bias function, to compute the variance of C \ Eq. (23) is approximated by 
var(C') = var model (C') + var noise (C" ) , (24) 

where var mode i(C') is defined by 

var model (C') = ] var model (C' , x)dx ) , (25) 

•'•min 


and 


var model(C",^) 
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0; otherwise 


(26) 


With C’ obtained by Eq. (18) and its variance given by Eq. (24), for a particular measurement at x, if the model 
is inadequate, we can state that, with about 68.2% confidence, the true value C + b(x) resides in the region of 

[c'-^/var(C'), C'+7var(C') . 



As an example to illustrate our theoretical results developed in this section, we fit our model Eq. (17) to a 
measured data set of v = 0.01 x x 2 + Norm(l) in the jt region of [- 25, 25] , where 0.01 x x 2 is the true value and 
Norm( 1) is the measurement error which is Gaussian distributed with a Gaussian width of one. We use both the 
equal-weight and the Maximum Likelihood approaches with the number of data points of 51 and 201 respectively. 

In Figs. 4c and 4d, we show our fitting results, along with the fitting results from the equal-weight approach 
shown in Figs. 4a and 4b. In the figures, the diamonds represent te measurement data, the dashed lines for the real 
values of the data, and the solid lines for the fitting results with the bars indicating the standard deviations of the fit- 
obtained C’. Figs. 4a and 4b show that the standard deviations obtained from the equal-weight approach are 
unrealistically too small, whereas the standard deviations in Figs. 4c and 4d are about correct demonstrated by 
C'±std(C') covering slightly more than half of the real values. 
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Fig. 4. Fits of Eq. (17) to O.Olx + Norm{ 1) over the x range of [-25, 25], using the equal-weight and Maximum 
Likelihood approaches, respectively with the number of data points of N=5 1 and 1V=201 . (a-b) Equal-weight approach, 
(c-d) Maximum Likelihood approach with model error handling. 


For a future study, we would like to extend our simple model Eq. (17) to more complicated ones and compute the 
effects of model errors to the coefficients in those models. It is also desirable to compute std(b ) more accurately. 


5. SUMMARY 



In this proceeding, we showed a Maximum Likelihood approach to determine the c-ratios in the quadratic 
polynomial detector spectral radiance model. We gave mathematical expressions to compute the weights used in the 
approach, considering the random noises from both the dependent and independent variables. In addition, we 
developed a methodology to compute the impact of model errors on the c-ratios. We applied our Maximum 
Likelihood method to the VIIRS RSB measurement data. Our results revealed that the obtained c-ratios are in good 
agreement with those computed with the equal-weight method even though the quadratic model is may not be 
inadequate over the radiance range specified by the requirement. Critical to the knowledge of the measured radiance 
error, our Maximum Likelihood with model error handling approach yielded much more realistic c-ratio variances 
that are much larger than those from the data regression generated covariance matrix. We pointed out that because 
the biases are not random numbers, the uncertainties of the fitting parameters can not be reduced simply by fitting 
more data points. Therefore, in the presence of significant model errors, the covariance matrix values from data 
regression routines incorrectly address the variances of the fitting parameters. The Maximum Likelihood approach 
shown in this proceeding may lead to better sensor calibration, including better understanding of sensor calibration 
uncertainty. 
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APPENDIX 


In this appendix, in order to evaluate our results given in this proceeding, we shall carry out a computer 
simulation to generate dn jn and dn out . We then and fit the simulated data to extract the c-ratios and their 
variances, following the procedure described in Section 2. The simulation assumes that the true radiance, after 
divided by q , follows a cubic polynomial in dir. 

^ = ^ + dn + ^xdn 2 +^xdn i . (A.l) 

C 1 c \ c l c \ 


In the simulation, we set co/q=-0.85, c 2 /ci=-3x10 6 , c 3 /q=- 1.6x10 9 ,and t = 0.566 . The results are 
shown in Fig. A1 and Fig. A2. From Fig. Al, we can see that the Maximum Likelihood approach fit-obtained r 
(shown by the solid horizontal line) is much different from the true r . Fig. A2(a) shows the results from applying 
the Maximum Likelihood method with model error handling. Fig. A2(b) shows the results from the equal-weight 
approach. The large difference between the retrieved and the true radiances is indicative of flaws in the attenuation 
method used in this study. 

To resolve the large difference between the retrieved and the true radiances, shown in Fig. A2, we carry out fits 
with the attenuator’s transmittance fixed at its true value of 0.566. The results are shown in Figs. A3 and A4. From 
Fig. A4, we can see that in this case the retrieved radiances are much closer to the real ones, indicating that we may 
need to know the attenuator’s transmittance before using the attenuation method. 
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Fig.Al. Results from the fits to simulated data, following the procedure in Section 2. Diamonds are calculated T and 
the horizontal solid line is for the fit-obtained r from the Maximum Likelihood approach. 
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Fig.A2. The ratio of the retrieved radiance over the true radiance minus one. (a) Maximum Likelihood approach, 
(b) Equal-weight approach. 
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Fig.A3. Results from the fits to simulated data, with t fixed at 0.566. Diamonds are calculated T and the 
horizontal solid line is for the fit-obtained x from the Maximum Likelihood approach. 
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Fig.A4. The ratio of the retrieved radiance over the true radiance minus one at r = 0.566 . (a) Maximum 
Likelihood approach, (b) Equal-weight approach. 
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